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Quantum-coherent coupling of a mechanical oscillator 
to an optical cavity mode 
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Quantum control of engineered mechanical oscillators can be achieved by coupling the oscil- 
lator to an auxiliary degree of freedom, provided that the coherent rate of energy exchange 
exceeds the decoherence rate of each of the two sub-systems. We achieve such quantum- 
coherent coupling between the mechanical and optical modes of a micro-optomechanical 
system. Simultaneously, the mechanical oscillator is cooled to an average occupancy of 
n=1.7±0.1 motional quanta. Pulsed optical excitation reveals the exchange of energy be- 
tween the optical light field and the micromechanical oscillator in the time domain at the level 
of less than one quantum on average. These results provide a route towards the realization 
of efficient quantum interfaces between mechanical oscillators and optical fields. 
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Mechanical oscillators are at the heart of many 
precision experiments, such as single spin detec- 
tion [1] or atomic force microscopy and can ex- 
hibit exceptionally low dissipation. The possibility 
to control the quantum states of such engineered 
micro- or nanomechanical oscillators, similar to 
the control achieved over the motion of trapped 
ions [5], has been a subject of longstanding inter- 
est [H H] , with prospects of quantum state trans- 
fer [SHS] , entanglement of mechanical oscillators [9] 
and testing of quantum theory in macroscopic sys- 
tems [101 111) . However, such experiments require 
coupling the mechanical oscillator to an auxiliary 
system — whose quantum state can be controlled 
and measured — with a coherent coupling rate that 
exceeds the decoherence rate of each of the sub- 
systems. Equivalent control of atoms has been 
achieved in the context of cavity Quantum Elec- 
trodynamics (cQED [T^) and has over the past 
decades been extended to various other systems 
such as superconducting circuits [12], solid state 
emitters [14] or the light field itself [15] . 

Recently, elementary quantum control at the 
single-phonon level has been demonstrated for the 
first time, by coupling a piezo-electrical dilatation 
oscillator to a superconducting qubit [16j. An al- 
ternative and highly versatile route is to use the 
radiation-pressure-induced coupling of optical and 
mechanical degrees of freedom, inherent to opti- 



cal microresonators |17j . which can be engineered 
in numerous forms at the micro- or nanoscale [18]- 
'20|. This coupling can be described by the inter- 
action Hamiltonian H = hgQa)a{w -\- b), where a 
(6) is the photon (phonon) annihilation operator, 
h is the reduced Planck constant and go is the vac- 
uum optomechanical coupling rate. In the resolved 
sideband regime (where the mechanical resonance 
frequency ftm exceeds the cavity energy decay rate 
k), with an intense laser tuned close to the lower 
optomechanical sideband, one obtains in the rotat- 
ing wave approximation the effective Hamiltonian 



H = hg{aP +a^b 
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(1) 



for the operators a and b now displaced by their 
steady state values. We have introduced here the 
field-enhanced coupling rate [HJ |^ g = y/n^go, 
where Uc denotes the average number of photons in 
the cavity. In the absence of decoherence, the uni- 
tary evolution ([T]) corresponds to swapping of the 
(displaced) optical and mechanical quantum states 
with a period of 27r/f2c, where fJc = Sg is the co- 
herent energy exchange rate. This state swapping 
is at the heart of most quantum control protocols 
[3"^. In practice, however, this unitary evolution 
is compromised by the coupling of both degrees of 
freedom to their respective environments. Hence, 
it is important for f^c to exceed both the optical 
decoherence rate k, and the mechanical decoher- 
ence rate 7, defined here as the inverse time needed 
for a single excitation to be lost into the environ- 
ment. Reaching this regime of quantum-coherent 
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FIG. 1: Optomechanical resonators. (a) 

Scanning electron micrograph of a spoke-anchored 
toroidal resonator with 31 ^m diameter used 
for the optomechanical experiments reported in 
this work (scale bar 10 /xm). (b) Simulated dis- 
placement (exaggerated for clarity) of the funda- 
mental radial breathing mode of the structure, 
(c) Equivalent optomechanical Fabry-Perot cavity: 
Quantum-coherent coupling is achieved when the 
enhanced coupling rate ilc is comparable to or ex- 
ceeds the optical and mechanical decoherence rates 
(k, Fmnni)- Due to the large asymmetry between 
mechanical and optical frequencies, the occupa- 
tions of the two baths are widely different. 



coupling has proven challenging. While the onset 
of normal mode splitting [5TJ [22] (which occurs for 
12c > '«/2, when k, largely exceeds the mechanical 
dissipation rate T^) could be observed in a room 
temperature experiment |23| . the mechanical de- 



coherence rate 7 = Fj, 



1 



(where 



n^ = k^T/hVt^) is strongly enhanced by the cou- 
pling to the environment at equilibrium tempera- 
ture T. The condition Jlc > (k, 7) is analogous to 
the strong coupling regime encountered in atomic 
cavity QED [T^] with the atomic decay rate be- 
ing replaced by 7, and the Rabi frequency by the 
optomechanical coupling rate fJc- Recently, a su- 
perconducting microwave cavity electromechanical 
system has accessed this regime [53J [5S]. If im- 
plemented at optical frequencies, the availability 
of well-developed quantum optical techniques, and 
the possibility to transport quantum information 
through room-temperature fiber links could give 
access to a range of new applications [5HH1 ^M ■ 

Here we demonstrate for the first time the 
quantum-coherent coupling of an optical cavity 



field to a micromechanical oscillator. The exper- 
imental setting is a micro-optomechanical system 
in the form of a spoke-supported toroidal optical 
microcavity [26] . Such devices exhibit high quality 
factor whispering gallery mode resonances (with a 
typical cavity decay rate K/27r < 10 MHz) coupled 
to mechanical radial breathing modes via radiation 
pressure. The vacuum optomechanical coupling 
rate qq = ^xzvm can be increased by reducing 
the radius R of the cavity. However, the larger per 
photon force Hto/R is then usually partially com- 
pensated by the increase in the mechanical reso- 
nance frequency rim — and correspondingly smaller 
zero point motion xzpm = \/^7(2woff^^in)- More- 
over, small structures also generally feature larger 
dissipation through clamping losses. To compen- 
sate these opposing effects we use an optimized 
spoke anchor design (cf. Fig. IT] and appendix) 
that maintains low clamping losses and a moderate 
mechanical resonance frequency while reducing the 
dimensions of the structure. Devices fabricated in 
this manner (with i? = 15 /im) exhibited coupling 
rates as high as go = 27r x 3.4 kHz for a resonance 
frequency of 78 MHz, as determined independently 
at room temperature [27] . 

To reduce the mechanical decoherence rate 7 = 
Fm^im, the microcavity, coupled to a tapered op- 
tical fiber, is embedded in a helium-3 cryostat 
|28| . For low optical power, the ■^He buffer gas en- 
ables thcrmalization of the resonator over the en- 
tire cryogenic temperature range (Tmin = 650 mK) 
in spite of its weak thermal anchoring to the 
substrate. A quantum-limited continuous-wave 
Ti:Sapphire laser provides both the coupling field 
and homodyne local oscillator used to measure 
the weak phase fluctuations imprinted on the field 
emerging from the cryostat by mechanical dis- 
placement fiuctuations. While the coherent cou- 
pling rate V,c can be determined unambiguously by 
probing the coherent response of the system [25], 
the mechanical decoherence rate is affected in a 
non-trivial way by the light-absorption-dependent 
sample temperature and the mechanical mode's 
coupling to its environment |28| . which is domi- 
nated by two-level fluctuators at cryogenic temper- 
atures [30]. In order to systematically assess the 
decoherence rate, the coupling laser's frequency 
wi = Wc -l- A is swept in the vicinity of the lower 
mechanical sideband. This allows to bring the dis- 
placed cavity mode a (of frequency |A|) in and 
out of resonance with the mechanical mode b (of 
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FIG. 2: Optomechanical interaction in the 
weak coupling regime (i^c < k/2). (a) The co- 
herent response of the system is obtained by sweep- 
ing a weak probe beam (created via phase modula- 
tion of the coupling laser at frequency rimod) over 
the cavity resonance and recording the homodyne 
signal. The fitted detuning is indicated for each 
of the traces. The fit provides accurate estimation 
of the coupling rate Qc via the Optomechanically 
Induced Transparency window (see appendix for 
details), (b) The measured Brownian noise spec- 
trum in the absence of a coherent probe for each of 
the dctunings in (a). The spectra are corrected for 
the detector response and a small contribution of 
Guided Acoustic Wave Brillouin scattering in the 
fibers (see appendix). The inset shows a close-up 
of the spectrum obtained for A/fim — —1.89. The 
indicated decoherence rates 7 are deduced from the 
amplitude of each noise spectrum, (c) Inferred me- 
chanical displacement spectra calculated using the 
extracted parameters, (d) Retrieved occupancy as 
a function of detuning. 



frequency O^). For each detuning point, we ac- 
quire the coherent response of the system to an 
optical excitation in a first step using a balanced 
homodyne detection scheme. These spectra (Fig 
2k) allow to determine all parameters of the model 
characterizing the optomechanical interaction (ap- 
pendix). For large detunings |A| > il^, they es- 
sentially feature a Lorentzian response of width 



K and center frequency |A|. The sharp dip at 
fl w i7,n originates from Optomechanically In- 
duced Transparency [29j, and for fim = —A, its 
width is approximately ^1/k. The coupling rate, 
as derived from a fit of the coherent response for 
a laser power of 0.56 mW, is fl^ = 2^/n^go = 
2tt X (3.7 ± 0.05) MHz (corresponding to an intra- 
cavity photon number of fie = 3 • 10^). 

Additionally, for each value of the detuning the 
noise spectrum of the homodyne signal is recorded 
in the absence of any external excitation (Fig. [2b). 
The observed peak represents the phase fluctua- 
tions imprinted on the transmitted light by the 
mechanical mode's thermal motion. The constant 
noise background on these spectra is the shot-noise 
level for the (constant) laser power used through- 
out the laser sweep (see appendix for details). 
The amplitude of the peak is determined by the 
coupling to and the temperature of the environ- 
ment, and therefore allows to extract the mechan- 
ical decoherence rate. All parameters now being 
fixed, it is moreover possible to retrieve the me- 
chanical displacement spectrum (Fig. [2];). As can 
be seen, for detunings close to the sideband, when 
the (displaced) optical and mechanical modes are 
degenerate, the fluctuations are strongly reduced. 
This effect of optomechanical cooling can be under- 
stood in a simple picture: In the regime flc ^ k, 
the optical decay is faster than the swapping be- 
tween the vacuum in the displaced optical field and 
the thermal state in the mechanical oscillator. In 
this case, the mechanical oscillator is coupled to 
an effective optical bath at near-zero thermal oc- 
cupancy finiin with the rate Fcooi = ^c/'^- Ideally, 
n-min = ^^/ISf^m ^ 1 is govcmed by non-resonant 
Stokes terms a^b^+ab neglected in the Hamiltonian 
(fTj). Any excess noise in the coupling beam will, 
however, cause an effective increase of finiin and 
precludes any quantum state manipulation due to 
the impurity of the field's quantum state. In prac- 
tice it has proven crucial to eliminate phase noise 
originating from Guided Acoustic Wave Brillouin 
Scattering [31 in the optical fibers by engineering 
their acoustic modes using HF etching (appendix). 

Evaluating the mechanical decoherence rate for 
A = — r^m at a cryostat setpoint of 0.65 K, we find 
7 = 27rx(2.2±0.2) MHz— significantly smaller 
than Qc- Simultaneously, the average occupancy 
of the mechanical mode is reduced to n = 1.7±0.1, 
which is limited by the onset of normal mode split- 
ting. Indeed, as f2c approaches k, the thermal 
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FIG. 3: Quantum-coherent coupling. The coherent optical response (a) and the incoherent me- 
chanical noise spectrum (b) for various detunings of the coupling laser (with constant power P = 1.4 
mW). All curves are normalized, and vertically displaced by the detuning. Evident in both panels is the 
avoided crossing which originates from optomechanical normal mode splitting. A second much smaller 
mechanical mode at 71 MHz is omitted, (c) Homodyne noise spectrum obtained for A = — ^im. The 
red line is a fit of the model. Only the decoherence rate (determining the amplitude) is fitted, while 
the shape is fixed through the parameters determined from the coherent response, (d) Comparison of 
the mechanical (green) and optical (blue) decoherence rates with the coherent coupling rate (red) as a 
function of detuning. The increase of 7 close to resonance reflects heating of the cavity due to the larger 
amount of absorbed light. On the lower mechanical sideband, where the interaction is resonant, the 
decoherence rates are comparable to the coupling rate, achieving quantum-coherent coupling. 



fluctuations are only partially dissipated into the 
optical bath, and partially written back onto the 
mechanics after one Rabi cycle. Note that this 
occupancy, corresponding to 37% ground state oc- 
cupation, is associated with a sideband asymmetry 
of 1 — ri/(ri -I- 1) « 40%. A measurement of this 
asymmetry, similar to those performed on trapped 
ions [2], would yield direct evidence of the quan- 
tum nature of macroscopic mechanical oscillators. 
We subsequently increase the strength of the 
coupling fleld to reach Vie ~ k. The signature of 
normal mode splitting pTVE^ can be seen from 
both the coherent response and the fluctuation 
spectra in Fig. [31 Both detuning series exhibit 
a clear anti-crossing, the splitting frequency be- 
ing 5.7 MHz. The decoherence rate (see Fig. [3]i) 
is slightly raised compared to Fig. |2] due to laser 
heating and a higher buffer gas temperature of 
0.8K, amounting to 7 = 27r x (5.6 ± 0.9) MHz at 
the lower mechanical sideband. We hence demon- 
strate r2c/7 = 1-0, which constitutes a four-orders 
of magnitude improvement over previous work in 
the optical domain |23j , and brings the system into 



the regime of quantum-coherent coupling. 

As a proof of principle, we flnally demonstrate 
the dynamical exchange of weak coherent excita- 
tions between the optical and mechanical degrees 
of freedom in the time domain (Fig. 111). To this 
end, we modulate the coupling laser's phase at the 
mechanical resonance frequency, with a Gaussian 
envelope of 54 nano-second duration. This modu- 
lation creates a pair of sidebands — one of which is 
resonant with the optical cavity — that contain on 
average ten quanta per pulse (cf. appendix). By 
measuring the homodyne signal it is possible to 
directly observe the coherent exchange of energy. 

In the regime of weak coupling (Fig. [4k) the 
optical pulse excites the mechanical mode to a 
finite oscillation amplitude, which decays slowly 
with the nearly unperturbed mechanical dissipa- 
tion rate. This in turn only imprints a weak 
signature on the homodyne signal. Increasing 
the laser power, we reach J7c/27r — 11.4MHz > 
K/27r = 7.1MHz, and the envelope of the ho- 
modyne signal — corresponding to the amplitude 
of the measured sideband field — undergoes sev- 
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FIG. 4: Coherent exchange between the 
optical field and the micromechanical os- 
cillator probed in the time domain as measured 
( "Data" ) and calculated numerically ( "Model" ) . A 
modulation pulse (blue traces) applied to the phase 
modulator creates a light pulse probing the dynam- 
ics of the optomechanical system. The response 
of the system is encoded into the optical pulse at 
the output of the coupling fiber as recorded by 
the homodyne receiver (red traces, 250,000 aver- 
ages). Using the full model of the system, the me- 
chanical displacement can be simulated in addition 
(green traces) . In the regime of weak coupling (top 
panel), the optical output pulse exhibits only a 
weak signature of the nearly unperturbed mechan- 
ical ringdown excited by the short burst of radia- 
tion pressure. In the case of strong coupling (mid- 
dle and bottom panel) the modulated envelopes of 
the time-domain response indicate several cycles of 
oscillation between (very-low energy n = 0.9, bot- 
tom) coherent optical and mechanical excitations 
of the system. 



eral cycles of energy exchange with the mechanical 
oscillator, before it decays with a modified rate 
of (k + Fni) /2 « k/2, corresponding to the de- 
cay rate of the optomechanical polariton excited. 
Using our model matched to previously taken co- 



herent response measurements, we can derive not 
only the homodyne signature — which reproduces 
our data very well — but also the expected mechan- 
ical oscillations resulting from this pulsed excita- 
tion (see Fig. l4| . These reveal very clearly how the 
excitation cycles continuously between the optical 
and mechanical modes. Although the employed 
detection is far from optimized for time-domain 
experiments, we have achieved a signal-to-noise ra- 
tio of 40 by averaging 250,000 traces within a to- 
tal acquisition time of approximately two minutes. 
Reducing the excitation further so that one pulse 
contains less than one photon on average, a clear 
signature of several swapping cycles can still be ob- 
served (Fig.jiJ;). Replacing the weak coherent state 
by a single photon [7 , we expect Rabi-like oscilla- 
tion of the Fock state from optics to mechanics and 
back. A repeated quadrature measurement yield- 
ing a bimodal distribution would be an unambigu- 
ous signature of the quantum nature of the state 
after the full swap |32j . 

Obtaining quantum- coherent coupling flc > 
(7, k) has several interesting consequences. First, 
it allows the mapping of in principle any quantum 
state of the optical field onto the mechanical mode 
via the use of a time-dependent coupling field. As 
a simple example, initialization of the mode in 
the ground state can be efficiently achieved in this 
regime using a 7r-pulse which swaps the thermal 
state of the oscillator and the vacuum in the opti- 
cal field [33] . Note that the manipulation of large 
quantum states becomes increasingly challenging 
since the lifetime of the number state \n) scales 
with 1/n. In this context, it will be beneficial 
to further reduce the decoherence rate by limit- 
ing spurious laser heating and employing materi- 
als with low intrinsic loss, as well as to increase the 
optomechanical coupling rate by further miniatur- 
ization. The regime of quantum-coherent coupling 
demonstrated here has been proposed as a gen- 
eral quantum link between electromagnetic fields 
of vastly different frequencies, e.g., different wave- 
lengths in the optical spectrum or microwave and 
optical photons '^,'M'. In that respect, the efficient 
coupling of the demonstrated system to a low-loss 
single mode optical fiber is beneficial. Moreover, 
quantum-coherent coupling enables the use of the 
mechanical oscillator as a transducer to link oth- 
erwise incompatible elements in hybrid quantum 
systems, such as solid-state spin, charge, or super- 
conducting qubits and propagating optical fields 
"SI. In this context, we note that the mechanical 



modes of silica microresonators have already been 
shown to be intrinsically coupled to effective two- 
level systems in the form of structural defect states 
in the glass [50] , 

In conclusion, the reported experi- 
ments — demonstrating quantum coherent coupling 
between a micromechanical oscillator and an op- 
tical mode — represent a first step into the 
experimental investigation and quantum optical 
control of the most tangible harmonic oscillator: 
a mechanical vibration. 
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Al. EXPERIMENTAL DETAILS 
A. Experimental setup 

Figure |A1| shows a schematic of the employed experimental setup. At the heart of the optical setup 
is a Ti:Sapphire laser (Sirah Matisse TX) operating at a wavelength around 780 nm. The laser exhibits 
quantum limited amplitude and phase noise at Fourier frequencies relevant for this experiment. During 
the experiments the laser is locked to an external reference cavity such that drifts of the laser detuning 
A can be neglected during the acquisition time. 

The sample itself resides in a helium-3 exchange gas cryostat (Oxford Instruments Heliox TL) that is 
used for cryogenic pre-cooling of the mechanical mode to low temperatures. Since the toroids are situated 
directly above the surface of the liquified helium-3, the achievable temperature is directly linked to the 
vapor pressure curve of helium-3. As the toroidal microstructures are thermally very well isolated from 
the substrate, one relies on cooling via the helium-3 exchange gas. As a consequence, cryostat temperature 
setpoints of at least 650 mK (corresponding to pressures larger than w 0.15 mbar) are favorable. Coupling 
of light into the toroid is achieved via a tapered optical fiber that is approached using piezo positioners, 
which are compatible with low temperature operation (Attocube GmbH). The fiber ends are guided 
through and out of the cryostat and constitute one arm of the interferometer that is part of the balanced 
homodyne detection scheme. The length of one of the ca. 8 m long arms is servo-locked using a movable 
mirror to cancel the DC-component of the interferometer's signal. This setting allows a shot-noise limited 
read-out of the phase noise imprinted onto the transmitted laser field. The laser passes an electro-optical 
modulator that allows to create sidebands around the laser frequency. Last, the laser power is stabilized 
actively in absolute terms at the input of the experiment to ensure operation at a constant light intensity. 



Ti:Sapphire 
laser 



Phase 
control 




FIG. Al: Setup. See text for details. 



The three colored building blocks highlighted in Figure XT depict the three different measurements 
that are routinely performed one after the other. 

• Coherent response. A network analyzer sweeps the upper modulation sideband over the optical 



resonance and demodulates the corresponding (coherent) signal (cf. section ASF I. 



• Noise spectrum. Connecting only an electronic spectrum analyzer gives access to the incoherent 
noise spectrum (cf. section A3E). 



• Time domain response. Sending a pulsed stimulus from an arbitrary waveform generator to the 
EOM which modulates the coupling laser gives access to the dynamic time domain response of the 
optomechanical system (cf. section AlC). 



B. Influence of guided acoustic wave Brillouin scatting 



A crucial prerequisite for optomechanical measurements in the quantum regime is the use of a quantum 
limited laser source. From the point of view of quantum manipulations, added noise in the coupling beam 
corresponds to an improper state preparation, the optical beam being in a statistical mixture of pure 
quantum states. In the weak coupling limit fie ^ k, where the optical field acts as an effective bath, 
these extra fluctuations correspond to an increased temperature of the bath and prevents cooling close 
to the quantum ground state [TH3]- In addition, classical laser noise driving the optomechanical system 
can lead to ambiguous signatures such as squashing in the noise spectra, as reported previously [1]. We 
have verified in our previous work [5j, that the employed laser source is quantum limited. However, as 
is well known from fiber-based quantum optics experiments |6j, optical fibers can give rise to classical 
phase noise, in the form of guided acoustic wave Brillouin scattering (GAWBS). This process involves 
thermally driven radial mechanical modes of the fiber, that also modulate the optical path length. 
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To investigate the presence of GAWBS we have recorded the noise spectrum from the homodyne 
detector when the fiber is retracted away from the cavity in an imbalanced Mach-Zehndcr interferometer. 
Several classical peaks are observed on top of the shot-noise background (not corrected for the detector 
response) (Fig. A2). The width (i.e. damping) of the noise peaks was observed to narrow dramatically 



when the buffer was partly stripped off the fiber, clearly demonstrating the mechanical nature of the 
peaks. One of the peaks coincides with the mechanical resonance frequency of 78 MHz. However, the 
frequency of the dilatational fiber modes is proportional to the inverse fiber radius and can therefore 
be shifted by etching the fiber cladding in an HF solution. Immersing the fibers (without removing the 
acrylate buffer, which is permeable to HF) in a 40 % HF solution for 50 minutes reduced the cladding 
diameter from 125 /im to 95 /im. This increased the GAWBS mode frequencies of all fibers in the setup 
by « 30%, shifting them away from the mechanical resonance frequency of the toroid (Fig. A2i). 

The data shown in Fig. |A2[ 3 have been taken under the same conditions as the lowest occupation run 
shown in Fig. 2 of the main manuscript. The residual noise at 78 MHz (due to small portions of fibers 
that have not been etched) is approximately a factor seven smaller than the initial peaks, corresponding 
to a noise level of approximately 2 % of the shot-noise. This noise is generated along the fibers, both 
before and after the cavity. Special care was taken to minimize the length of unetched fiber before the 
cavity. The fact that an influence of the cavity detuning and coupling parameters on the transduction of 
these classical noise peaks into a measured signal is not discernible indicates that indeed most remaining 
noise originates from fiber after the cavity. Under this assumption, the independent noise of the GAWBS 
can be subtracted from the signal in order to estimate the decoherence rate and occupation. Figure 
A2 ; in the main text shows that the shape of the spectra, as predicted from independently measured 
parameters, is in excellent agreement with the data after subtraction, in which no signs of squashing are 
observed. Nonetheless, we have performed an additional analysis for the lowest-occupancy data under 
the assumption that half of the noise is generated before the cavity, which leads to deviations of the 




b) 



40 60 80 100 120 140 
Frequency [MHz] 



1.0 



i 




4' 




— before 

— after etching 

— shot noise level 


^^ 









70 



75 80 

Frequency [MHz] 



85 



90 



FIG. A2: Engineering of the fiber GAWBS noise spectrum. The figure shows a broadband 
background spectrum of the imbalanced homodyne signal where the GAWBS modes are visible. Blue, 
green and red trace are taken with unmodified, partly stripped buffer and (almost entirely) etched fiber 
respectively. As expected for guided dilatational acoustic waves of the optical fiber the frequencies are 
increased by a factor of about 1.3 for a thinned fiber of around 95 /im diameter (as compared to 125 /im 
before). Doublets in the red trace are due to slightly different final etching radii (difference is about 
3;um) of the different fibers in our setup (i.e., local oscillator fiber and the signal fiber). The difference 
in relative heights of the peaks is attributed to varying readout conditions, and as such only the peak's 
frequencies are of interest. The inset shows a zoom of the background for the final setup (i.e., etching 
reduced diameter fibers, balanced homodyne arm lengths) and illustrates the achieved improvements, 
i.e., the reduced contribution of GAWBS to the background at the mechanical resonance frequency. 



11 

decoherence rate and occupation of 7% and 5%, respectively. This upper bound of the influence of 
GAWBS, corresponding to an uncertainty of 0.08 phonons, is included in the quoted errors (cf. section 



A3II 



C. Time-domain response 

In order to probe the coherent dynamics of the optomechanical system in the time domain, the strong 
pump beam is tuned to the red sideband, and an RF pulse, resonant with the mechanical oscillator, is 
sent to the EOM. The upper modulation sideband excites the strongly coupled system. The subsequent 
evolution of the transmitted signal is recorded using the homodyne detector and an oscilloscope. 

An arbitrary signal generator (Agilent 33250A) is used to generate the RF pulses. The time dependent 
voltage U(t) is a sine wave modulated by a Gaussian envelope: 

U{t) = E{t) sin(f],„odi + M (Al) 

E{t) = C/oe-(^)' (A2) 

with a carrier frequency ilmod — 2-k x 77 MHz and an envelope duration r = 32 ns (FWHM = 54 ns). 
A digital oscilloscope, synchronously triggered with the signal generator is used to record and average 
the homodyne response. The very small signal originating from the balanced detectors is amplified and 
filtered, around a frequency of 75 MHz, with a bandwidth of 100 MHz. For low excitation amplitude, 
averaging is necessary to extract the coherent response out of the incoherent thermal and quantum noises 
from the optomechanical system. 

The modulation depth j3{t) corresponding to the instantaneous value of the slowly varying envelope is 
given by: 

/3(0=^^, (A3) 

where Vt^ = 154 V is the voltage corresponding to a tt phase shift of the beam in the EOM (NewFocus 
4002). For a weak modulation depth (/3 ^ 1), a fraction (/3/2) of the optical carrier power Pc is scattered 
into each of the two first modulation sidebands. The total optical power in the upper sideband can hence 
be simply approximated by 

m.P.(.fl)' (A4) 

The total energy in the pulse can then be obtained by integrating the instantaneous power over the 
duration of the pulse. The average number of photons in one pulse is hence given by: 

A2. OPTIMIZED SPOKE-ANCHORED TOROIDAL RESONATOR 

A. Sample design 

The optomechanical microresonators investigated in this work are specially designed toroidal micro- 
cavities, optimized to achieve large optomechanical coupling rates and small dissipation. Toroidal silica 
whispering gallery mode microresonators exhibit mechanical modes coupled to the optical modes through 
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radiation pressure [7_. Of particular interest is the lowest order radial breathing mode (RBM), whose 
motion maximally modulates the optical cavity length. In the context of quantum-coherent coupling, 
it is important to simultaneously achieve large values of ^dl ^-^d ^c/n, where Qc = 2go |a|- The 
vacuum optomechanical coupling rate go is given by ^ \/ti/ {^mcsflm): where R is the toroid radius, 
rUeff is the effective mass, and lu and fini are the optical and mechanical resonance frequencies, respec- 
tively. For a given incident power, optical frequency, and environment temperature, and assuming the 
resolved sideband regime and the coupling laser being tuned to the lower mechanical sideband, one ob- 
tains ric/7 oc ^R^^^y^ltn^rles/K] and Hc/k cx I i?ilin y^mcifK ) . It is therefore obviously beneficial 

to reduce the sample dimensions to decrease R and rriee- However, such miniaturization is generally ac- 
companied by an increase of the mechanical frequency fl^, as well as an increase of T^ due to enhanced 
clamping losses. In microtoroids, both of these adverse effects can be countered in a design in which the 



toroid is suspended by spokes from the central pillar, as shown in Fig. |A3^ . 

The introduction of spokes serves three purposes. First, they isolate the mechanical motion of the 
toroidal RBM from the pillar support, strongly reducing clamping losses [S]- Second, they reduce the 
mechanical mode volume and thereby the effective mass. Third, the effective spring constant is reduced, 
which lowers the mechanical resonance frequency. In practice, one needs to carefully consider the precise 
dimensions and positioning of the spokes, as these strongly affect both clamping losses Fdamp and go. 
Figure [XSJ d shows the displacement profile of the RBM of a spoke-supported toroid of radius i? = 15 /Ltm 
for various combinations of spoke length and position, as simulated with a finite element method. The 
Si02 thickness is 1 /im, the minor toroid radius is 2 fiia, the spoke width is 500 nm, the pillar diameter 
is 1 /im, and the toroid is vertically offset from the middle Si02 disk by 400 nm. Since we are interested 
in the RBM only, it suffices to simulate 1/8 portion of the microresonator while assuming symmetric 
boundary conditions on both of the two 'cut' planes. As can be seen from these examples, the mechanical 
mode profiles can change drastically depending on the spoke dimensions. Of the examples in Fig. \A3}p , 
only 'D' depicts a mode that is purely localized to the outer toroid, with purely radial displacement, as 
illustrated in the cross-section in Fig. |A3[:. 



The origin of this wildly varying nature of the RBM is revealed in Fig 
the inner disk (defining the spoke placement) and the spoke length Is are varied systematically. 
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FIG. A3; Sample optimization by finite element modeling. See text for details. 
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colorscale depicts the parameter F, defined as 



2ttE, 



mcch 



F^ —■"""' _ (A6) 

cpniJ^JAz{^)\'dA 

Here, i?moch is the total mechanical energy in the mode, c the speed of sound in silica, p the density 
of silica, and Az (x) the out-of-plane displacement amplitude, with the integration extending over the 
area Ap of the interface between the pillar and the silica disk. F is proportional to the expected value 
'-'^ ^ciamp' 'when the clamping area Ap is considered as a membrane radiating energy with a power 
P = cpfl'^^J^ |A2;(x)| dA |8j. A previous study has experimentally found a correspondence of F « 

(SFciamp/ (2vr))^ for larger toroids. As can be seen from the figure, the expected clamping losses vary 
strongly with spoke dimensions, ranging from 10"'^ to 10^ Hz. Most notably, several lines can be identified 
in this parameter space where clamping losses are large (indicated by the dashed lines). For parameter 
combinations along each of these lines, the REM frequency approaches that of another mechanical mode 
of the structure. As a result, the two modes exhibit an anticrossing, with the hybridized modes showing 
a character of both uncoupled modes. This is the case for examples 'A', 'B', and 'C in Fig. |A3| 3, which 
show the RBM hybridized with a fiexural mode of the inner Si02 disk, the outermost Si02 membrane, 
and the spoke itself, respectively. In the vicinity of these anticrossings, the vertical displacement at the 
pillar, and as such the radiation into the substrate F~^, are strongly enhanced. To achieve a design that 
exhibits small clamping losses, it is therefore crucial to avoid these parameter regions, as is the case for 
mode 'D' in Fig. [XSj j. 

The aforementioned anticrossings affect the coupling rate go as well, albeit to a lesser degree. Fig. 



A3 3 shows go, calculated as in [9], assuming the optical mode is localized at the edge of the toroid with 
negligible transverse size. At the anticrossings, go is reduced (i.e., the effective mass is increased), as a 
significant part of the mode's energy is in that case associated with displacements that do not modulate 
the cavity length. Away from the anticrossings, however, the RBM mode is localized exclusively in the 
toroid and outermost part of the membrane, well isolated from the inner disk and pillar support. As a 
result, rricB is nearly identical to the physical mass of this volume. It is therefore important to minimize 
the volume of the outermost membrane, i.e., the distance between the spokes and the toroid. As can 



be seen in Fig. XS', this simultaneously allows to reach the smallest possible resonance frequency. In 



practice, the laser reflow process used to form the toroid poses a lower limit on the remaining distance 
between spokes and toroid. 



B. Sample fabrication 

To fabricate the spoke-anchored microresonators, we use a combination of optical lithography and dry 



etching techniques outlined in Fig. A4 In a first step (b), a disk including the spokes is transferred in a 
1 /im thick film of thermal oxide on a Si wafer (a), through optical lithography followed by reactive ion 
etching of the Si02. In a second photolithography step (c), smaller disks of photoresist are defined that 
cover the center of the Si02 disks, including the spokes. These serve to protect the exposed Si surface 
between the spokes during the subsequent isotropic XeF2 etch (d) of the Si substrate. Care is taken 
to stop the etch shortly before it reaches the apertures in the Si02 disk. After removing the protective 
photoresist disks, a laser refiow of the underetched disk is performed (e), forming the silica toroid. Finally 
(f), a second XeF2 etch releases the toroid and reduces the pillar diameter, typically to a value smaller 
than 1 fim. 
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^^ photoresist i^bi Si02 ^^ Si 
FIG. A4: Sample fabrication. See text for details. 

C. Sample characterization 



The vacuum optomechanical coupling rate qq is measured at room temperature in a vacuum chamber. 
Therefore, the mechanical motion is read out using an external cavity tunable diode laser at 1550 nm, 
that is locked to a cavity resonance. In order to avoid any radiation pressure effects we perform these 
measurements at very low laser power (typically around 100 nW). The transmitted light is amplified by 
a low noise erbium-doped fiber amplifier and sent onto a photodector. For absolute calibration of the 
mechanical spectrum registered by an electronic spectrum analyzer, we use a phase-modulation technique 
[TO] . We extract a vacuum optomechanical coupling rate of g^ = 1700Hz for a wavelength of 1550 nm 
(i.e. go = 3400Hz at 780 nm). 

The mechanical linewidth measured at room temperature (8.1 kHz) was found to be higher than ex- 
pected from the calculated F-parameter (cf. section A2 A). However, performing the same measurements 
in the cryostat (cf. Fig. ESp), a linewidth as low as 3.6 kHz was found on the same microresonator, indi- 
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FIG. A5: Sample characterization, a) shows a calibrated mechanical noise spectrum for the resonator 
used throughout the manuscript (except the last panel of Fig. 4 of the main manuscript), which was 
measured at room temperature in vacuum. The fit (red line) was used to extract a vacuum optomechanical 
coupling rate of 3400 Hz at 780 nm. b) mechanical damping of this toroid vs. cryostat temperature. The 
red line is a fit according to the TLS model presented in [51 [TT] , the grey lines represent the contributions 
from resonant (dotted) and relaxational (dashed) processes. The fit yields a negligible contribution of 
clamping losses. 
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eating that a loss mechanism other than elamping losses must dominate at room temperature. Sinee there 
losses due to two level fluctuators (TLS, [TTJ[T2]) have been found to be significantly lower (Hnewidths 
below 4 kHz have been measured at room temperature for conventional toroids of similar frequency) , we 
believe that the dominating loss mechanism is thermo-elastic damping (TED) (THj. At low temperatures, 
where TED is strongly reduced, the main loss mechanism is coupling to TLS. Figure |A5[ 3 shows the 
measured temperature dependence of the mechanical linewidth at low temperature, obtained with the 
laser (with 100 nW power) resonant with the, in this case, strongly overcoupled optical resonance to 
avoid dynamical backaction. The variation of T„^ with temperature can be fitted using a model for the 
TLS losses [5J [TT] . It is found that this mechanism dominates the total losses for all reachable cryogenic 
temperatures. This means that it is not possible to retrieve an accurate estimation of the temperature- 
independent contribution Fdamp- We can however conclude that it must be at least smaller than 2 kHz 
for this sample. This shows that in our optimized spoke-supported design, we have successfully mitigated 
the clamping losses to the level where they are insignificant compared to intrinsic dissipation. 

A3. MODELING OF OPTOMECHANICAL INTERACTION 

This section summarizes the theoretical model which was used to extract all entities of interest from 
our data. Figure [A6| shows the parameters and variables of the model, and their mutual connections. 



A. Conservative dynamics 

The conservative dynamics of an optomechanical system are described by the Hamiltonian |14) 

H = -hn^ [q^ + p^) + hwc ia^ a + -\ + hgoq a^a, (A7) 

where mechanical quadrature operators q and p are related to the corresponding mechanical ladder 
operators b and b^ via 

q = b + b'f (A8) 

p^{b~b^)/i. (A9) 

With these definitions [q,p] = 2j, and the actual mechanical displacement and momentum are given by 
x' = XzpfQ and p' — hp/2xzpF with the amplitude of the zero-point motion 



^ZPF = J ^ ^„ . (AlO) 

The vacuum optomechanical coupling rate go quantifies the strength of the optomechanical interaction 
and is given by go = ^o^^zpf with go = dujc/dx. 

B. Quantum Langevin equations 

The Hamiltonian ( |A7[ ) determines the conservative evolution of the optomechanical degrees of freedom. 
Optical and mechanical dissipation, and the corresponding fiuctuations, can be taken into account by 
introducing the mechanical dissipation rate Fm and the optical dissipation rate n = kq -\- Kcx (where 
Kex represents losses to the coupling waveguide and kq all other optical losses) as well as the optical 
noise terms ^sin, Sscnv and the thermal Langevin force, which we express as a rate Sfth by writing the 
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FIG. A6: Theoretical model used. See text for details. 

physical force in raomentum units of /i/2xzpf. This leads to the well-known Langevin equations of cavity 
optomechanics P^HTT] 

d{t) = (lA - -J a{t) - igoq{t)a{t) + ^/^^(sin + 6sin{t)) + ^/Ko<5scav(0 (AH) 

q{t) = rt^p{t) (A12) 

m = ^nnMt) - 2g^a\t)a{t) - T,^p{t) + <5/th(t), (A13) 

where the convention A — lj\ — lj^ was used to denote the detuning of the laser (angular) frequency u)\ 
from the bare cavity resonance frequency ujc, and a is expressed in a frame rotating at uj\. 

In order to accurately model the response of the optomechanical system over a wide range of parameters 
(detuning, Fourier frequency, optical and mechanical excitation) for a single set of parameters, we have 
refined this generic model by including other effects which arc known to be inherent to most optical 
microcavities, and are discussed in the following. 

Photothermoelastic backaction. Thermoelastic forces driven by temperature gradients induced by light 
absorption can induce mechanical displacements. The starting point to model these displacements are 
the coupled equations of motion known from the standard theory of thermoelasticity |13j 



^ V^u + (A + A*) V(V • u) + / = (3A + 2/i)aV6l + pii 

fctV^S — ctpO — (3A + 2/i)aro(V ■ u) — v KabsO' a. 



(AM) 
(A15) 



These equations connect the displacement field u{r, t) and the temperature elevation 0(r, t) above the 
mean temperature Tq. Here, A and p are the Lame parameters, a the thermal expansion coefficient, p 
the mass density, / a body force (e.g. due to radiation pressure), kt the thermal conductivity, Kabs the 
photon absorption rate, ct the heat capacity, and the function v{r) describes the spatial distribution of 
the light absorption. Evidently, a thermoelastic body force 



/to(r,t) = -(3A + 2M)aV0(f,t) 



(A16) 
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acts on the mechanical modes of the structure when a temperature gradient \J9{f,t) is present. In the 
scalar representation of the mechanical dynamics, we therefore have to add a thermoelastic force ftcit) 
proportional to the material parameters A, ji and a, as well as an overlap integral of the mechanical 
mode's displacement pattern and the temperature gradient V9{f^t). Assuming that the temperature 
gradients are predominantly driven by the absorption of laser light in the resonator, one can express the 
scalar photothermoelastic force as 

/ptc(0 = Xptc(i) * K.ahsa^{t)a{t), (A17) 

where we have absorbed the spatial overlap integrals between the mechanical and (the gradient of) 
the thermal modes, as well as the thermal modes and the spatial pattern of light absorption into the 
magnitude of the function Xpte(i)- The temporal dynamics of the adjustment of the relevant temperature 
gradients to a changing amount of light absorption is represented by the time-dependence of Xpte(i) ("*" 
denotes a convolution). Note that while this formulation accounts for the quantum fluctuations of the 
intracavity field a(i), the statistical nature of photon absorption events is neglected. This is justified 
considering that the quantum fluctuations of optical heat deposition ( "photothermoelastic shot noise" ) 
have a much smaller effect on the mechanical mode than the direct fluctuations of the radiation pressure 
term 2go a^a. 

Dynamic photothermorefractive frequency shift. A temperature elevation 9{f,t) within the optical 
mode volume furthermore changes the refractive index, and therefore the optical resonance frequency. In 
analogy to the description in the previous section, we are using a simple scalar description of the form 

Aa;ptr(i) = Xptr(i) * K^^hsO^ {t)a{t) (A18) 

for this frequency shift, where the response function Xptr{t) accommodates spatial overlap integrals of 
the light absorption pattern v{r) and the thermal modes as well as the temporal dynamics of the latter, 
and, in addition, the spatial sampling of the induced refractive index changes 

An{f,t) ^ ^9{r,t) (A19) 

by the optical mode. 

Thermorefractive noise. The local temperature elevation 9{f,t) also undergoes thermal fluctua- 
tions — independent of the presence of light. Within a volume V, they amount to squared fluctuations of 

ra 

(«K.f>, = ^. (A20, 

The spatial distribution of these fluctuations can be calculated using a Langevin ansatz [19 , by adding a 



fluctuational source term to the heat diffusion equation (A15). Predominantly via the thermorefractive 
effect (dn/dT ^ 0) , the resulting temperature fluctuations again induce resonance frequency fluctuations 
SuJti-{t). Its temporal correlation function (or equivalently, power spectral density) have been estimated 
for simple whispering-gallery mode resonator geometries 



Taking these additional three effects into account, the equations of motion can be written as 
a{t) = (i(A - Aujpt.it) - Sujtrit)) - ^) a{t) - igoq{t)a{t)+ 

+ VK^iSin + Ssin{t)) + y/K^SSci,vit) (A21) 

q{t) = n,Mt) (A22) 

p{t) = ~n,^q{t) ~ 2goa^{t)a{t) - r„,p(t) + SftiM + fptc{t). (A23) 
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C. Linearized model 

A large coherent field sent to the optomechanical system induces a relatively large classical intracavity 
field a, and induces a displacement of the mechanical mode by q. If the system is stable around this 
steady-state, the dynamics of the small fluctuations around this equilibrium are described by a set of 
equations obtained via the substitution a{t) = a + Sa{t) and q(t) = q + Sq{t), and retaining only first-order 
terms in the fluctuations. This yields 

Sa{t) — ( +iA — — j Sa{t) — iKabsaXptr(^) * {a*Sa{t) + aSa' {t)) — igoa6q{t) — 

- iaduJt,.{t) + y/K^SSin{t) + y/K^SScav{t) (A24) 



n-^^ [dq{t) + r„, Sm + nl Sq{t)] = -2goidda^it) + d*6ait)) + 

+ Sftbit) + XptcW * {dSa^it) + d*Sa{t)) (A25) 

with A = cji — (ujc + goq + Kabs|ap(Xptr(0 * !))■ This set of equations is best solved in the Fourier 
domain, yielding 

(— i(A + rj) + k/2) 5a{il) = — iKabsflXptr(^)(fl*^fl(^) + dSa* (il)) -~ igQd5q{n)+ 

— idSu!tr{fl) + -y/Kox'5sin(^) + \/KO'^*cav (^) (A26) 

- ^-^5q{n) - (-2.90 + Xptc(f^)) {d6aHn) + d*6a(n)) + SUin). (A27) 



For simplicity, we refer to the Fourier transform of the respective functions by simply writing them with 
a frequency {fl) argument. Note that 6a^{il) denotes the Fourier transform of 5a^{t), equal to [(5a(— il)]^^; 
and that [(5(7(— £7)]^^ — 5q{il) for the Hermitian operator Sq{t). 

To further simplify the problem, we approximate the response functions of the photothermal effects by a 
single-pole, low-pass response, assuming implicitly that the relevant temperature (gradient) distributions 
adjust themselves only with a certain delay to a change in the absorbed optical power. Assuming that 
this delay is larger than the relevant oscillation periods considered here, one can approximate 

Xpt.(fi)«^^ (A28) 

Xptc(fi) « 25ptc^ (A29) 

and finally obtains 

{~i{A + n) + k/2) 6a{n) = dgptr^{d*Sa{n) + dSa^il)) - igQdSq{n)+ 

— idSuJti{fl) + y/UcK^Sini^) + y/KoSSca.vi^) (A30) 

-02 - tnT„, +J^^^(f^) = _2 (go + ^9pte^) (dda^n) + d*6a{n)) + Sft^{n). (A31) 



On 



These equations are used to calculate the coherent response and fluctuation spectra (cf. sections A3E 



ASF) 



D. Homodyne detection 

The optomechanical experiment is embedded into one arm of a balanced homodyne interferometer. 
At the initial beamsplitter, the laser field (and fluctuations in the fiber mode) are split up into a 'local 
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oscillator' arm, and the arm that serves as input to the cavity: 



Sin = \/l - rsias - Vrshs (A32) 

sio = \A-sias + Vl - rshs, (A33) 

evidently valid both in time and frequency domain. Here, we also take into account the vacuum fluctua- 
tions Sshs entering the beamsplitter at the unoccupied port, 

Slas = Slas + ^Slas (A34) 

Sbs = ^Sbs- (A35) 

The field Sjn drives both the mean field a and the field fiuctuations within the cavity, as described in the 
previous section. The intracavity field a, in turn, couples back into the single-mode fiber taper, and the 
usual input-output formalism gives the field Sout at the output of the cavity via the relation 

Sin Sout == V '*oxO (AoDJ 

We furthermore take into account that only a fraction T^cryo of the light power at the output of the cavity 
is measured as 'signal' in the homodyne detector due to optical losses, e.g. in the cryostat. For yycryo < 1, 
we again have to account for quantum vacuum ^Scryo that enters the optical mode, 

Ssig Y^ciyoSout ' \ ^ ^cryoScryo V '^ ' / 

Scryo ^ ^Scryo- (^Aooj 

Finally, in the homodyne receiver, the differential signal 

6h = s-ioe+'^'°,54g + si;e-^*'°<5s,ig + s^i^e-'^^'Ssl + sl.^e+'^^'Ssi, (A39) 

is measured. The fiuctuational terms 6h and 5q of interest can then be expressed as a linear function of 
the fiuctuations driving the system, 

[Sq] =^^'('^*las (^sfas ^^hs Ssl^ 5Scav <5sJav (^Scryo ^sj^y^ SuJtr Sfth) ■ (A40) 



Here, the coefficients of the matrix M follow directly from the relations (A30l-( A39) 



E. Calculation of noise covariances 



We assume that all input noise terms of eq. ( A40[) can be described by zero- mean Gaussian noise 



operators whose variances are known. Representing the covariances between two noise operators x and 
2/ as a symmetrized spectrum Sxy{^) defined according to 

i {{x{n), y{n')}) = 2ns,y{n) 5{n + n'), (A41) 

the only non-zero covariances are characterized by the spectra 

for the optical quantum noise entering the system, 

Ssf,,sUn) ^ 4n^T^ (A43) 
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for the thermal Langevin force, where we have assumed ri-ni ~ k^T /Ml-a^ 3> 1, and 

S5^,,Su,A^)=St,nm, (A44) 

for the thermorefractive noise |20j . whose contribution we found to be neghgible in the data presented 



in this manuscript. By the Unearity of equation (A40), it fohows that the covariance matrix iVout of the 



output noise operators is then related to the input covariance matrix TVjn by the simple expression 

^out = M{+n) ■ TVin • M{-^f. (A45) 

F. Coherent dynamics of the system 

In order to calculate the coherent response of the system to the probing by a phase-modulated input, 



eq. ( A40 ) can be used. By assuming a sufficiently narrow detection bandwidth and/or sufficiently large 



phase modulation of depth 5lp, one can set 

(SsJav ~ ^Scav ~ 5sl^ sa 6shs ~ 5sl^yo ~ (^Scryo ~ ^/th ~ ^^tr ~ (A46) 

5sus = isus^V, (A47) 

and calculate the frequency-dependent transfer function from a phase modulation 5^ to the homodyne 
signal 5h. 

This coherent response is obviously directly measured in the sideband sweeps that we routinely perform 



(cf. section Al A|. Moreover, by multiplication of the (complex) spectrum of the excitation pulse with this 
transfer function, the response of the homodyne signal in the time domain can be numerically determined 
via the inverse Fourier transform. 



G. Analysis of the coherent response 

The coherent response spectra are important to accurately extract the different parameters of the 
optomechanical interaction as well as to calibrate the mechanical noise spectra. A typical coherent 
response is shown in Fig. |A7^ . The Lorentzian peak centered around Omod — 140 MHz results from 
the absorption of the upper modulation sideband by the cavity and reflects the optical response of the 
system. The maximum of the homodyne signal is obtained when the modulation sideband is resonant 
with the cavity. Hence, the center frequency and width of this peak correspond to the detuning |A| 
and the linewidth k, of the cavity, respectively. The sharp feature at rimod = ^m is the manifestation of 
Optomechanically Induced Transparency [23| ; an interference effect due to the resonant excitation of the 
mechanical mode. For weak coupling power and/or large detuning, the dynamics of the mechanical mode 
is hardly affected by the optomechanical interaction and the width of the dispersive feature is given by 
the mechanical linewidth T^. 

For larger laser power, the width of the OMIT window increases, reflecting the width of the damped 
mode Fni-f ri^K/ (k^ -|- 4(A -t- f^m)'^) • Hence, the fit of the coherent response allows to extract the coupling 

rate Vt^ and the corresponding intracavity field a. We introduce ag = a/ _.^\_ ,„ to obtain a parameter 



-iA-l-K/2 



independent of detuning. The model of eq. (A40), assuming pure radiation pressure backaction, fits the 



measurements well (cf. Fig. A7i). However, as can be seen in Fig. A7d, a small systematic deviation 
appears for high coupling power. This systematic eflfect is very well reproduced by the model including 
the photothermoelastic effect. 

Finally, Figure |A7| : shows a series of coherent response spectra taken for decreasing laser detunings, 
and a laser power of 0.6 mW. The observed increase of the amplitudes for small detuning can be fitted 
accurately by introducing the photothermorefractive effect in the model (red lines). The parameter 
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FIG. A7: Fitting the model to the coherent response, (a) A coherent response spectrum taken 
with a power of 0.56 mW, at T=0.65 K. (b) Spectrum for 1.4 mW, at T=0.8 K with fits including the 
photothermoelastic effect (red line) and without (yellow dashed), (c) Spectra for 0.6 mW at T=0.75 K, 
for various detunings. The photothermorefractive effect is included in the fitted model and accounts for 
the increased amplitude for small detuning. 



^ptr introduced here is dependent on detuning since the thermorefractive coefficient ^ depends on 
temperature. 

For example, we have extracted from the fits to the full detuning series of Fig. 2 of the main manuscript 
f7,n/27r = 78.2MHz, K/27r = 6.0MHz, oq = 14.2 • 10^ (with go/2'K = 3.4kHz), 5ptc/27r = -122Hz and 
gpti./27r = 0.32 Hz (at the lower mechanical sideband). 



H. Extraction of the decoherence rate 



The fits of the coherent response spectra determine all parameters characterizing the optomechanical 
interaction, and therefore the transduction of mechanical displacement fluctuations to optical fluctuations. 
The spectral shape of the noise originating from the Langevin force is thus fixed, so that the amplitude of 
this contribution can be fitted using the model of eq. ( A45 1 . As wc fit the spectral density of the actually 



measured voltage signal, these extracted amplitudes depend on the gain of the subsequent detection chain, 
which is not precisely known. 

This ambiguity is removed by a calibration technique [Sj based on a reference phase modulation, which 
allows to relate noise spectra taken under arbitrarily different acquisition conditions. In this manner, 
we link the low-temperature noise spectra to a measurement at a higher cryostat temperature (4 K), in 
which a high helium gas pressure, and low optical power (~ 100 nW) ensure the thermalization of the 
sample, so that the Langevin force is known to an estimated accuracy of 3%. In this high-temperature 
measurement, a known phase modulation is applied, whose amplitude can be compared with the coherent 
response spectra acquired with every low-temperature measurement. Assuming that no drift occurs in the 
phase modulation chain, this method allows to absolutely calibrate the Langevin force — and therefore 
the mechanical decoherence rate — in the low-temperature measurements. Importantly, this derivation 
reveals possible changes of the decoherence rate both due to a changed temperature (bath occupation 
TT-ni) and mechanical dissipation rate T^. 



22 

I. Error analysis 

We use the large number of traces acquired during a detuning sweep to estimate an error on each 
of the four parameters assumed to be independent of the detuning (flm, k, ao, ffpte)- This is achieved 
by successively letting each of these parameters vary with the detuning, while the three others are 
still fitted globally. The error on each parameter X is obtained by calculating the standard deviation 
AX = \/{{X — -^^0)^)1 where Xq is the value obtained when all four parameters are kept constant over 
the whole detuning range. Advantageously, this procedure reflects also systematic errors due to drifts 
of the experimental settings over the detuning series, and physical effects that are not captured by the 
model. The following uncertainties were obtained with this method for the run presented in Fig. 2 of the 
main manuscript: an/27r = (78.2260 ± 0.0007) MHz, K/2n = (6.04 ± 0.08) MHz, Sq = (14.2 ± 0.2) x 10^, 
(7ptc/27r — (122 ±52) Hz. These errors, by affecting the shape of the expected noise spectra, also translate 
in an error on the fitted decoherence rate and occupation. A Monte-Carlo approach is used to assess the 
final error on 7 and h. The fit of the noise spectrum is repeated with a set of randomly drawn parameters, 
assuming an independent normal distribution for each of the previous parameters. Importantly, the 
resulting uncertainty depends on the particular detuning point. On the lower optomechanical sideband, 
the standard deviation of the results is given by ( A7/7)jjjqJj,[ = 6 % and (^"-/^)modGi = 4 %. Another 
source of uncertainty for these two parameters is the independent calibration of the optomechanical 
transduction that we estimate to be on the order of Acaiib = 3% from the scatter between calibration 
measurements taken at different probing power. Finally, as discussed in section |A1 B[ an uncertainty 
Agawbs is quadratically added to account for the possible presence of GAWBS in the optical fibers 
before the cavity. The total error for this example is given by 




Acaiib + AgA^gg — 10 % 



odd 



Acaiib + A^AWBS - ''"^«- 
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